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Abstract 

We construct exact solutions for a system of two nonlinear partial differen- 
tial equations describing the spatio-temporal dynamics of a predator-prey 
system where the prey per capita growth rate is subject to the Allee effect. 
Using the (^) expansion method, we derive exact solutions to this model for 
two different wave speeds. For each wave velocity we report three different 
forms of solutions. We also discuss the biological relevance of the solutions 
obtained. 
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1. Introduction 



In this paper we derive exact solutions for a diffusive predator-prey system 

Ut = u xx — f3u + (1 + f3)u 2 — u 3 — uv 

vt = v xx + kuv — mv — 5v 3 (1) 

where k, 5, m and are positive parameters. Subscripts x and t denote 
partial derivatives. The equations are expressed in dimensionless variables, 
where scaling of space and time have been introduced so to have the equation 
in a simple form. The biological meaning of the each term has been discussed 
in which we briefly review. The model is of the predator-prey kind, with 
u and v being the densities of prey and predator. Spatial redistribution of the 
population is governed by diffusive dynamics, with both species having the 
same diffusivities, j^]. In absence of the predator, the temporal dynamics of 



the prey density is of the Allee type 



0,3, 



small populations being not viable. 



Presence of the predator population negatively affects the prey population. 
In turn, the predator population is totally dependent on prey availability, 
as the only term in the predator equation that represents the growth of the 
population is the kuv term, where k quantifies the gain in natality due to prey 
consumption. The parameters and m represent the per capita mortality 
rate of prey and predator respectively, in the linear, small populations, limit. 
Finally, the 5v 3 is a closure relation taking into account the effects of higher 



trophic levels 



a. 



To investigate the dynamics of the above diffusive predator-prey system 



the authors of Ref. [l] have assumed the following relations between the 



parameters, namely m = (3 and k + 4= = ft + 1. In other words it has been 
assumed that the per capita mortality rate of prey and predator are equal 
and the rate of biomass production at the lower level must be consistent 
with the rate of biomass assimilation at the upper level of the web js, 9, 10 ]. 
Under this assumption Eq. (JT]) reads 

ut = u xx — flu + (k H — i=)u 2 — m 3 — uv 

v 5 

Vt = v xx + kuv — (3v — 5v 3 . (2) 

In our work also we consider the predator-prey system (T5]) only. We construct 
exact analytic solutions to the equation (j2J) in order to understand the prop- 
erties of this model for different parametric values. To do so we employ the 



(J^r) expansion method ll|, [l2|, |_13|, [lj, |l5|. This method has been applied 
to several nonlinear evolutionary equations. Here we demonstrate that the 
utility of this method in exploring dynamics of a diffusive predator-prey sys- 
tem. While implementing (^jf) expansion method to equation (J2J) we obtain 
exact solutions for two different wave speeds (vide Eqs. (I14p and (fT5|) ). In 
both the cases we give three different types of exact solutions. 

The plan of the paper is as follows. To begin with, in Sec. 2, we describe 
the (^) expansion method. In Sec. 3, we consider Eq. §2$) and derive exact 
solutions of it. Finally, we present our conclusion in Sec. 4. 



2. The (^)-expansion method 

In this section we discuss briefly the method of finding exact solutions 
for a system of nonlinear partial differential equations (PDEs) using (^|-)- 
expansion method. 
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Suppose that the system of nonlinear PDEs is of the form 

P(u, v, u t , u x , v x , u u , u xt , u xx , ...) = 0, 
Q(u,v,v t ,u X} v X} Utt,u xt ,u XX} ...) = 0, (3) 

where u = u(x,t) and v = v(x,t) are two unknown functions and P and Q 
are polynomials in u = u(x,t) and v = v(x,t) and their partial derivatives. 
The (^t) expansion method involves the following four steps. 
Step 1 : Let us introduce the travelling wave reduction 

u(x, t) = «(£), v(x, t) = u(f ), £ = x-cb, (4) 

in the PDE Q so that the latter becomes 

P(u, v, —cu', u', v', c 2 u" , — cu", u", ...) = 0, 

Q(u, v, —cv', u\ v', c 2 u" , —cu" , w", ...) = 0, (5) 

where prime denotes differentiation with respect to the new variable £. 

Step 2 : Suppose that the solution of (jSJ) can be expressed by a polyno- 
mial in (^r), that is 

( G'\m /G'\rn-1 , Q ' \™-2 

u{£) = a rn[-g) +a m ^-) +a m _ 2 ^J +..., 
/G'\ n /G'\ n - 1 /G'\ n ~ 2 

«(£) = (6) 

with G = is the solution of the second order linear damped harmonic 
oscillator equation 

G" + AG" + fiG = 0. (7) 

In the above a m , (3 n , a m _i, a , A), A, and /1 are constants and a m 7^ 0, 

(3 n 0. The positive integers m and n can be determined by substituting 



in (jSJ) and considering the homogeneous balance between the highest order 
derivative and nonlinear terms appearing in (jSJ). 

Step 3 : Substituting fl6]) in (jSJ) and eliminating the variable G" in the 
resultant equations by using ([7]) one gets two polynomials equations in (^). 
Now equating each coefficients of (^r) m and (^j)™ to zero one obtains a 
set of algebraic equations for the parameters a m , p n , ee m -i, Pn-i---, a o, Po, A, 
and fi. Solving these algebraic equations one can get exact values for these 
coefficients. 

Step 4 : Substituting the values of a m , (3 n , a m _i, /3 n _i..., ao, /3o, A, and 
and c and the general solution of (J7J) in ([6]) one can obtain three different 
types of travelling wave solutions for the given system of nonlinear PDEs. 

3. Diffusive predator-prey system 

In this section, we apply the method described in the previous section to 
the nonlinear PDEs and construct exact solutions. Substituting into 
(T2J) we get the following system of ordinary differential equations (ODEs), 
namely 

u" + cu — Bu + (k H i=)u 2 — u 3 — uv = 0, 

v" + cv + kuv — (3v — 5v 3 = 0. (8) 

Suppose that the solution of ODEs (jSj) can be expressed by a polynomial in 
(^r) which is of the form ([6]). 

Substituting (jBJ) and their derivatives in OH]) and performing the homoge- 
neous balance between u" and u 3 and v" with in resultant equation we 
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find m = 1 and n = 1. So we fix the polynomials be of the form 

w(0 = ai(§-)+oo, v(0 = ft (§■)+#>, ai,/3i^0. (9) 

Substituting the expressions (J9J) and their derivatives in (IS]) and rearranging 
the resultant equation in the descending powers of we arrive at 

/ G' \ 3 ct 2 / G' \ 2 

[2cti - cKi] ^— J + [3«iA - cat + ka\ + — ^ - 3a 2 a - 

+[(2/i + A 2 )«i — cA«i — + H -p- — 3«q«i — «i/3o — «o/?i] 



a 2 

+(/i«iA — c/i«i — /3«o + «o "I — ?= — a o — a o0o) = 0, 



V6 



(10) 



[2/3i - ^i 3 ] (§) 3 + [30i A - eh + kaifr - 35f3 2 1 f3 ] 
+[(2fi + A 2 )/3i - cA0i - 00! + A;a o 0i + kaxfo - 35 

+(/i0iA - c/i0! - 00 o + A;a o 0o - <*0o) = 0. (11) 

Equating the coefficients of (^j) m , m = 0,1,2,3, to zero in equations (flOjl 
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and (11 ip we get the following set of algebraic equations, namely 



2cti - a\ = 



a? 



3aiA — c«i + ka\ H — -L — 3a 2 a — ai/3\ = 0, 
v 5 

(2/j, + A 2 )«i — cA«i — f3ai + 2A;«o«i H t= 3«q«i — ai/3 — a [3i = 0, 

V 5 

/i«iA — c/iai — /3a + «q H — -p — «q — a /3o — 0. (12) 

v 5 

2fr - = 0, 

3/3iA - c/3i + fcaiA - 3J/3f/3 = 0, 

(2/i + A 2 )/3i - cA/3i - (3/3i + ka ^ + fcai/3 - 35/3^/3i = 0, 

/i/3iA - c/i/3i - /3/3 + ka f3 - 5/3 3 = 0. (13) 

Solving the above system of algebraic equations (TT2"j) and (TT3"j) we obtain two 
sets of values for the constants ai, ao, A) and c: 

(a) ai = ±V2, (3 = — , Px = ±y-, c = =F-^=, A = =F — -/=— , 

/3 = A;a - «q + 2//. (14) 



c = ±i= (2fc - 3a + ^) , I = - H - 2fl){ ' k T + al ~ 2/i) ■ (15) 

Since both the sets separately satisfy the algebraic equations in ffT2l and ffTBl 
they individually form a compatible solution. From each set we derive an 
exact solution for the nonlinear PDEs (j2]). 

To begin with let us take the values given in ( IT4")) . With the values given 
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in ( EH]) ) the solution (JH]) reads 

«(« = ^g) +ao , „ (0 = ± ^g) + - , 16) 

It is known that the linear damped harmonic oscillator equation ([7]) admits 
three different types of solutions depending on the values of A and /i, namely 
Case 1: A 2 - 4/i > 

G(0 = U sinh ^ X2 ~ 4 ^ + C2 cos h ^FE±A {17) 

V 2 2 / 

Case 2: A 2 - 4/i< 

g(« = ( Cl C08 V^EZj + C2 8in ^EZ ? ) (18) 

Case 3: A 2 - 4/i = 

G(0 = (ci + c 2 0e ( - A/2)f • (19) 

Substituting ( TTT|) - (TT9|) into (Tl6|) we arrive at the following form of solu- 
tions, namely 

Case 1: A 2 - 4/i > 



_ ^ ( A | y^4^ / Cl cosh Vy^e + c 2 sinh 

V 2 2 U sinh ^£ + c 2 cosh ^ 

2/ A VA 2 - 4/i / Cl cosh ^Ef^ + C2 S inh 



w(0 = ± 



Of, 



*\ 2 2 v Cl sinh ^e+ C2 cosh ^f^yy V^' 

(20) 
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Figure 1: The densities of prey (solid line) and predator (dashed line) as given by the 
exact solution (|20|) shown for time t = 0; when A>2 V //I; Parameters are cto = 1.2, A; = 5.9, 
5 = 3,/^ = 0.2, ci = 20, c 2 = 10. 
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Figure 2: The densities of prey (solid line) and predator (dashed line) as given by the 
exact solution (|2T|) shown for time t = 50; when A<2y / /I; Parameters are ao — 3, k = 12.2, 
6 = 2, /x = 5, ci = 20, c 2 = -10. 
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Figure 3: The densities of prey (solid line) and predator (dashed line) as given by the exact 
solution (|22p shown for time t — 10 ; when A = Parameters are 6 = 3, k = 2.03, 

/i = 1, ci = 20, c 2 = 10. 

Case 2: A 2 - Afi < 



M(0 = ±V V 2 + "Sr" ( CTT ■ 1 ! + 

\ \ C\ cos — 4 + c 2 sin -^—2 — 4 

... ^ /2/ A , yi^ Z -csm^fg^ + ccos^fg^ , , 



«0 



(21) 



Case 3: A 2 - ifj, = 



where a = — 2v/2 ^ +fc ; ^ = x ± (^jVj c i an d c 2 are arbitrary constants. 
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For the second set of values we end up with the same form of solution 
with the only difference in the values in the parameters c, (3 and A: 
Case 1: A 2 - 4/i > 

, J A v^r^/ Cl cosh^^e + c 2 sinh^^A^ 

m= , [2( A y^^ / c.cosh^^e + c.sinh^^A ^ c 
V{) ~ V A"2 + 2 Usmh^ + c.cosh^J; ^ 

(23) 

Case 2: A 2 - 4// < 



M(0 = ±V ^I - 2 + 2~ I • ' ' ' 

\ ci cos g — 4 + c 2 sin g — 4 



, , _. A - A 2 / -ci sin + c 2 cos ^ A ^ 



1 Ci cos 2 4 + c 2 sin - y — g 4 ' 



(24) 



Case 3: A 2 - 4/i = 



where a = y/2jl, £ = x =F \ {^sfik — 3y/2a + ^f^j , C\ and c 2 are arbitrary 
constants. 



4. Discussion and Conclusion 

In this paper, we have constructed exact solutions for a diffusive predator- 
prey system which is modeled by a system of two coupled nonlinear PDEs. 
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Using the (^) expansion method, we have derived exact solutions for two 
different wave speeds. 

The solutions that we have obtained are singular and cannot be taken 
at face value describing actual situations in ecology. Notwithstanding they 
present a very interesting property: depending on the sign of A 2 — 4/i , the 
nature of the solution changes from a single structure to a periodic one, which 
is akin to pattern forming systems. 

We express A 2 — 4/x in terms of the original parameters in the equation, 
we get A 2 — 4/i = k/2 — 2(3. Therefore, if k 2 > 4/3 we have a single structure 
like in Fig.(l). And if k 2 < 4f3 we have a periodic structure. The meaning 
of k is that it measures the gain in natality obtained by the predator, and (3 
is its mortality. It follows that periodic structure formation comes from the 
strength of mortality. In a more intuitive way, /3 is the inverse of the typical 
time for a predator population to decay in absence of preys. If this time is 
short, we have a periodic pattern, the population continuously decaying and 
recovering. If this time is long, a smoother dynamics shows up. 

The equations for which we could find new solutions, besides the previ- 

I 

ously known, have the obvious drawback that matching of coefficients is 
necessary. This is the same situation in most cases in the subject of exact 
solutions for reaction-diffusion equations, beyond the specific cases of interest 



in biology, [16| . However, the broad view insight gained remains of interest as 
the systems considered contain many elements of more realistic, non solvable 
ones. 
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